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We study the stability and single-particle properties ol Fermi liquids in spatial dimensions 
greater than one via bosonization. For smooth non-singular Fermi liquid interactions we obtain 
Shankar's renormalization-group flows and reproduce well known results for quasi-particle lifetimes. 
We demonstrate by explicit calculation that spin-charge separation does not occur when the Fermi 
liquid interactions are regular. We also explore the relationship between quantized bosonic exci- 
tations and zero sound modes and present a concise derivation of both the spin and the charge 
collective mode equations. Finally we discuss some aspects of singular Fermi liquid interactions. 
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O ; I. INTRODUCTION 

o 

Landau's Fermi liquid theory is an early example of what we would now call bosonization. The anticommuting 
operators which appear in the bare Hamiltonian describing the interactions among fermions disappear in Landau's 
effective theory. Instead only c-number quasi-particle occupancies appear in the semi-classical energy functional. That 
the low energy semi-classical behavior of the Fermi liquid can be described in terms of these commuting variables 
suggests that a fully quantum bosonic description is obtainable. n 

Indeed, the Fermi liquid state itself is an example of a zero temperature quantum critical fixed pointtl. This fixed 
point is characterized by infinite U(l) symmetry which is not exhibited by the bare Hamiltonian. Thejnfinite U(1A 
symmetry simply reflects the conservation of quasiparticle occupancy at each point on the Fermi surfaceou. Shankarcl 
t— I ' has used the functional renormalization group (RG) approach to show that the Fermi liquid state is a generic feature 
of interacting fermions, at least at weak coupling and in the absence of the usual superconducting and charge- and 
spin-density wave instabilities. By establishing rigorous bounds, other workers have studied the stability question, at 
all orders in perturbation theory, but under more restrictive conditions such as a perfectly circular Fermi surface!! 
Jh ■ Haldane has asserted recently that a fully quantum description of Fermi liquids in dimensions greater than one is 
obtainable via bosonizationo. This viewpoint has been elaborated on by two of usa. In the present paper we continue 
, to develop this theory first by showing that Shankar's renormalization group result is obtained easily in the bosonized 
picture. Next we investigate the bosonic excitations in more detail. We show that collective modes are obtained 
in a semi-classical limit; furthermore the calculation of the single-particle boson Green's function yields information 
O ■ about the quasiparticle properties. In particular, by using the bosonization transformation to determine the fermion 
quasiparticle propagator we obtain well-known results for the fermion self energy: the imaginary part is proportional 
to cj 2 m|w| in two dimensions and just to 2 in three dimensions. We emphasize that the bosonization method yields 
non-perturbative information, so a natural next step would be to use it to study the effects of singular interactions. 
^ • We comment on the nature of two such singular interactions. 



II. RENORMALIZATION GROUP ANALYSIS IN THE BOSONIC BASIS 



Now that we know how to bosonize Fermi liquids we may use this picture to investigate the stability of the 
zero-temperature Fermi liquid fixed point to perturbations. First we reproduce the renormalization group results of 
Shankarcl in the bosonic basis. Three channels of fermion two-body interactions are marginal in the RG sense: forward 
scattering zero-sound (ZS), exchange scattering (ZS'), and Cooper pairing (BCS). For simplicity we consider a system 
of spinless fermions in two dimensions and a circular Fermi surface. The second assumption eliminates the possibility 
of nesting instabilities in the zero-sound channels which might produce charge or spin density waves. We also assume 
that the BCS coupling function Vbcs(S) is rotationally invariant. The BCS interaction pairs particles of equal but 
opposite momenta. For now we turn off the two zero sound channels; later we show that these channels have no effect 
on the renormalization of the BCSpinteractions. 

Fermi fields ip may be expressedtl in terms of the boson fields (f> as: 

x) = J= ^ e* k - x exp{z^(S; x)} O(S), (1) 
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where the dependence on time, t, is included implicitly in the spatial coordinates x. S runs from to 2n and labels 
the patch on the Fermi surface with momentum ks. V is the volume of the system which equals L 2 in two dimensions; 
the factor of V^ 1 ^ 2 is introduced to keep the fermion anticommutation relations canonical. Both the if) and <\> fields 
live inside a squat box centered on S with height A in the radial (energy) direction and width A along the Fermi 
circle. These two scales must be small in the following sense: kp >> A >> A. We satisfy these limits by setting 
A = kp/N and A = kp/N a where < a < 1 and N — ► oo. The quantity a in the bosonization formula Eq. ([!]) is 
a real-space cutoff given by a = 1/A. Here fl = A(JJ2tt) 2 equals the number of states in the squat box divided by 
A. Finally, O(S) is an ordering operator introduceddD to maintain Fermi statistics in the angular direction along the 
Fermi surface. (Anticommuting statistics are obeyed automatically in the radial direction.) 

With this connection between the fermion and boson fields we may check a number of relationships. For example, 
the fermion fields obey canonical anticommutation relations: 

{^(S;x) , V t (T;y)} = ^,T^(x-y) (2) 

because the boson fields in configuration space obey commutation relations: 

[0(5; x) , 0(T; y)] = ~ f! 2 6 s ,t e(h s ■ [x - y]) ; \x x - y ± \ « 1/A 

= 0; |&L-yj_|»l/A. (3) 

Here _L denotes directions perpendicular to the surface normal ris at patch S, and e(x) = 1 for x > 1; otherwise it 
equals -1. Normal ordered charge currents are defined in configuration space in terms of the both the Fermi and Bosc 
fields as: 

J(S; x) = V : ^(S;x)^(S;x) : 

= V lim {^(S; x + en s ) */>{S; x) - {^(S; x + en s ) ^(S; x))} 



= V4tt n s -V0(5;x) . (4) 
The momentum-space charge current is defined by: 

j(S; q) = J2 e ( s ; k + q) k ) {^I; q ^ - ^k} (5) 

k 

where 8(S; k) = 1 if k lies inside the squat box of dimensions A x A centered at S and equals zero otherwise. Given 
this definition, plus the fact that the Fermi fields in momentum and real space arc related in the usual way to preserve 
the canonical anticommutation relations with conventional normalization, 

t/>(S; x) = JL 0(S; k) e 4k ' x <&(k) , (6) 
v ^ k 

the two currents are related by a Fourier transform: 

J(S;x)=£ J(5;q) . (7) 
q 

Both curt cnts Eq. (§J) and Eq. (|) are dimcnsionlcss. The free Hamiltonian, written in terms of the Fermi fields, may 
also be bosonized using Eq. ([!]) and the result is quadratic in the <j> fields: 



H o = v F ^2j d 2 x ^ f (5;x) |^Z-jfe F | i/>(S;x) 

/rf 2 x{(n 5 .V)0(5;x)} 2 . (8) 
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The unusual prefactor of 4 ^ appearing in the bosonic Hamiltonian compensates for the anomalous right hand side 
of the boson commutation relations, Eq. (||), and thereby reproduces the correct spectrum. 

The next step is to bosonize the BCS interaction. To simplify the following algebra we set the Fermi velocity equal 
to one (vp = 1). A fermion in patch S of the Fermi surface is paired with a fermion in patch —S which is directly 
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opposite patch S. (Note that in angular coordinates, patches S and — S correspond to 9 and 9 + it, not 9 and —6.) 
The BCS action expressed in terms of the four Fermi fields is: 



Si„ t [ip, V^] 



E 

S,T 



dt d 2 x Vbcs( ; S T) tf(-S;x) ^t( 5;x ) ifj(r- x ) V(-r;x) . 



(9) 



Here, S and T only range over half of the Fermi surface to avoid double counting the pair interactions. The dimcn- 
sionless coupling function Vbcs must change sign under inversion because the fermions are spinless (Pauli exclusion 
principle) so Vbcs{6) = — Vbcs{^ + n )i a l so the interaction must be Hermitian so Vbcs(S — T) = Vbcs{T — S). To 
avoid sign mistakes it is important to keep track of the order of the fermion operators during the transformation to 
bosons. Formally the correct sign is set via the ordering operator O(S) but in practice it is easier to determine the 
sign by direct inspection of the fermion operators. To bosonize the interaction, each Fermi field is replaced by the 
right hand side of Eq. ([!]) which converts the interaction into the exponential of four <fi fields. 



'2vr' 



E 

S,T 



,,,2 Vbcs(S-T) 

at a x : ~ cos 

hp a z 
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(&r(x)-6r(x)) 



(10) 



where ^s(x) = 4>(S; x) + (f>(—S; x). 

Before implementing the RG transformation, first we discuss scaling at the zero loop level. Since we are concerned 
with scaling in the direction parallel to the surface normal rig, the integral over x space should be factorized into 
separate integrals over directions perpendicular and parallel to the Fermi surface normal. So d 2 x — dx± dx\\ with only 
dx» changing under scale transformations. Thus: A — * A, A — * A/s and a — > s a with s > 1. Clearly, dxu — > s dx\\ 
and dt — > s dt. The boson field 4> is invariant under the scale transformation as this leaves both the quadratic part of 
the action, So, and the BCS interaction invariant (marginal). 

Now we perform the renormalization group transformation on the BCS interaction to derive the flow equation. 
The fast parts of the 4>(T; k) fields are integrated out of the functional integral. To be precise, modes with momenta 
A/2s < |k • | < A/2 will be eliminated. In practice since it is easier to carry out the calculation in real space, instead 
we integrate out fields over short distance scales 2a < x • fiy < 2s a. Next, we rescale space and time: n^-x — > s fix-x 
and t — > s t. After performing these two operations we obtain the new BCS interaction coefficients Vbcs(S ~ T;s) 
and we may repeat the process. 

The integration over the fast modes is accomplished via the usual functional integral: 



exp(-S[0; a]) = J! / V ^ x ) **P(So ~ S*. 

t-> J 2a< |x-nr I <2sa 



(ii) 



Since only the fast modes are integrated out it is convenient to break the boson fields into two parts, the slow 
modes tf>' and the fast modes Now we may express the right hand side of Eq. ( pi] ) as (exp(— Si n t)) where 
the contraction is performed only over the fast fi> fields. The renormalized interaction is obtained by treating the 
interaction perturbatively and expanding (exp(— Si n t)) in powers of Vbcs- The first non-trivial term arises at second 
order, \{S 2 nt ), and since the interaction is diagonal in real-space, it is readily evaluated with the use of the real-space 
boson correlation functionO: 



(/ftTjx) mo)- f{T-M c £Lln( 



4-7T " X • flT + %T 

oo ; |xj_A| >> 1 



) ; \x^k\ « 1 



(12) 



Here and in Eq. ( |il| ) we have made a Wick rotation to imaginary time to avoid the poles along the real-time axis. 
The evaluation of the cosine-cosine correlation function that appears in \(S 2 nt ) is carried out by decomposing terms 
of the form 



47T , 



*>s(x)-<Mx)) 



into exponentials involving the slow and fast fields, and then using the identity: 

1 



e A e B 



= :e A+B : exp(AB + ^(A 2 + B 2 )) 
Since V BC s(S - T) = V BC s(T - S) we obtain: 



(13) 



(14) 
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(SL) = (^) 3 7T~~T2 I dT d%x E cos \^7T WrM - E V »csiR - T) V BCS (T S) 

77 ^ Fa ' J R,S L J T 



7T 



oo ./2a 



x — / dz^ / cZit|| — ^ o + {irrelevant operators} , (15) 



where uii is a spatial variable in the direction parallel to the surface normal at patch T and v is an imaginary time 
variable. The second step is to replace all the variables with the rescaled ones. This procedure does not change Eq. 
( |l5| ) since Vbcs(S — T) is marginal. Thus the f3- function is: 

dV BCS (R -S;s)_ a £ Vbcs{r t . s) Vbcs{t ^ s) (i6) 



d\n(s) V k F 



T 



The equation may be diagonalized by a Fourier transform over the interval [0, 2w): V m = J ^ e"™ e Vbos(^) ■ Note 
that only odd modes appear due to the requirement Vbcs(9) = ~ Vbcs(@ + 7r )- Then: 

^ =~\VI, (17) 



dln(s) 2 

which agrees with Shankar's result in the fermion basis at one-loop orderEi Clearly, a BCS instability exists if any of 
the channels are attractive (V m < 0). If all the channels are repulsive, the Fermi liquid fixed point is stable. 

We now turn on the other two marginal interactions, forward and exchange scattering, and ask whether the Fermi 
liquid interactions f c (S — T) which involve the ZS and ZS' channels alter the RG flows. 

H ^H + ±J2 I d2x fc(S-T) ^t(5 ;x ) V ,(s ;x ) Vj t(T;x) V(T;x) 

S,T J 

= H ° + ^ E / d " x M S - T ) [^-V0(5; x)] [n T -V0(T; x)] . (18) 

S,T J 

Unlike the BCS interaction, the forward and exchange interactions are quadratic in the boson fields and therefore 
parameterize different Gaussian fixed points, each with the infinite U(l) symmetry. This symmetry is reflected in 
the fact that the Hamiltonian Eq. (|l8| ) is invariant under changes in the phase of the fermions by different amounts 
in each patch: ip(S;x,t) — > e l9 ^ il>(S;x,t). Here we see the advantage of the bosonic representation: the Fermi 
liquid parameters are incorporated in a non-perturbative way into Hq. We carry out the calculation of the modified 
bosonic correlation function in the next section; here we just note that these modifications are sub-leading corrections 
to scaling that do not influence the leading RG flows of the BCS interaction. For example, though bosons in different 
patches are now correlated, this correlation is only of order times that of correlations within the same patch. So 
the leading behavior exhibited in Eq. (|l^) is unchanged. Thus we have the remarkable result that the leading-order 
stability of the Fermi liquid fixed point against Cooper pairing is unaffected by the existence of either small or large 
Fermi liquid parameters. n 

Actually, there is a sub-leading order instability: the Kohn-Luttinger effectcl. The bare V m due to, say, a short-range 
repulsive interaction are all positive but tend rapidly to zero at large-m. The ZS and ZS' channels, on the other hand, 
generate irrelevant contributions in the BCS channel (down by a positive power of A/kp) which renormalize the bare 
BCS interaction and therefore can make some of the V m slightly negative (unstable) at sufficiently large-m. Because 
of the small size of these coefficients, this effect is important only at extremely low temperatures and therefore we 
expect that the essential physics remains controlled by the Fermi liquid fixed point. 



III. INTERACTING BOSONS IN THE SEMICLASSICAL LIMIT 



Now that we have seen that the Fermi liquid fixed point is stable in the absence of attractive BCS interactions, 
we turn to the problem of diagonalizing the bosonic Hamiltonian that describes the fixed point. As we shall see, 
the problem is not as simple as it might seem at first because the current operators behave as both creation and 
annihilation operators. So we begin with an approximate semiclassical solution that bypasses this difficulty and yields 
the familiar collective mode equation. 

The particle-hole excitations of Fermi liquids have a bosonic character. Furthermore, the excitations are of either 
the charge or spin type: the bosonized Hamiltonian may be written as H = H c + H s to exhibit this factorization into 



4 



charge and spin sectors. The charge sector in D-dimensions (the volume V — L D now) is described by a Hamiltonian 
that is bilinear in the current operatorsa J(S; q): 



2 

S,T q 



H - = \ E E ^( s > T ; «i) J ( s ; -°0 J ( T ; <0 > ( 19 ) 

where fl = {^) D ^ 1 {^)- The Fermi liquid interactions are / C (S,T) = F C (S, T)/jV(0), where the density of states at 
the Fermi surface, summed over both spin species, is given by N(0) = for the case of the two dimensional Fermi 
gas. These interactions are incorporated into V c as matrix elements that couple currents in different patches: 

V C (S, T; q) = \ ft" V S^ 1 + I / C (S - T) . (20) 

Note that with this definition, and given the relationship between the currents and the <fi fields, Eq. the charge 
Hamiltonian H c of Eq. ( jig] ) agrees (up to a factor of 2 due to spin) with the form we found in the previous section, 
Eq. (|l8|). The charge currents obey the U(l) Kac-Moody relations: 

[J(S; q) , J(T; p)] = 2 S^ 1 6% +Pt0 Q q ■ n s ; (21) 

this algebra can be derived either from Eq. (||) and Eq. (||) or directly from Eq. (||) with the use of the canonical 
anticommutation relations for fermionsu. The quadratic form of this Hamiltonian implies immediately that it describes 
a fixed point invariant under the scale transformations A — > X/s. Similarly, the spin-sector is described by: 

H > = \ E E w T ; <o J ( s ' • J ( T ; ( 22 ) 

S,T q 

where the spin currents commute with the charge currents and obey the more complicated SU(2) Kac-Moody algebra!: 

[J Q (S; q) , J b (T; p)] = \ S^ 1 5 ab il q ■ n s <^ +p , + i ^J 1 e abc J C (S; q + p) (23) 

({a, 6, c} = {x,y,z} label the three components of the spin). Fermi liquid spin-spin interactions f s appear in the 
Hamiltonian as coefficients that couple spin currents in different patches: 

V S (S, T; q) = ~ v F (S) Q- 1 S^ 1 + i / S (S, T) . (24) 

Here we consider only the case of Fermi liquid interactions which are independent of the wavevector q. For regular 
interactions, q dependence only gives rise to additional irrelevant operators which do not change the behavior of 
the system at leading order. With singular interactions, on the other hand, divergences arise as q — » and these 
divergences may in some instances introduce relevant interactions that destroy Fermi liquid behavior. 

The equations of motion for the charge and spin currents yield the corresponding collective mode equations in the 
semi-classical limit. Using the Heisenberg equations of motion and the U(l) Kac-Moody algebra we readily obtain: 

^J(S;q) = [J(S;q), H c ] 

= V F q.n s J(S;q) + q.n s — £ / C (S - T) J(T; q) . (25) 

T 

for the charge sector. The first term on the right hand side has its origin in the free dispersion relation for particle-hole 
pairs of momentum q at patch S. The second term couples currents in different patches. Note that 

so the second term reduces to the usual integral over the Fermi surface in the N — > oo continuum limit. The equation 
of motion for the non-abelian spin currents contains, in addition to these two terms, a third term which makes the 
spins precess, when the system is magnetically polarized, in the local internal magnetic fielda: 
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z^J a (S;q) = [J a (S;q), H a ] 



v F q • n s J a (S; q) + \ q • n s - £ /,(S - T) J a (T; q) 



2 " ° y 

T 

- ^ E ^ k ) E ^ S - T) J C (T; q k) (27) 

k T 

Note that the factor of 2/3 multiplying the Fermi velocity vp in the free part of Eq. ( p4| ) does not appear in Eq. 
(p7|). The origin of the 2/3 factor is easy to understand in the Sugawara construction of the free fermion Hamilto- 
nian out of current bilinearsEI: it reflects the SU(2) invariance of the spin currents which permits the replacement 
J(S; q) • J(S; — q) — > 3J Z (S; q) J Z (S; q) for the purpose of computing the spectrum. The derivation given here, on the 
other hand, does not rely on this argument as spin rotational invariance is respected explicitly. Rather, the factor of 
2/3 in Eq. ( pi| ) cancels contributions to the free spectrum which arises from both the S ab and the i e abc terms in the 
non-Abelian anomaly. 

When an external magnetic field is applied, and |q| is small, the third term in Eq. (|27]) dominates. In the 
opposite limit of zero applied magnetic field, the spin equation is identical in form to the charge equation. Evidently 
bosonization captures all of the physics of charge and spin collective modes. The derivation is straightforward and 
relies only on the existence of the quadratic Hamiltonian and the Abelian and non-Abelian Kac-Moody commutation 
relations. In fact this approach may be useful in the study of highly spin-polarized Fermi liquids, a problem which 
has been examined recently by Meyerovich and Musaelin via the Green's function approachE3. 

As it stands these operator equations are exact, at least in the N — > oo limit in which the Kac-Moody algebras Eq. 



( |2l| ) and Eq. (23) become exact. The difficulty in obtaining exact solutions to either of these equations originates in 
the fact that neither algebra is equivalent to the canonical commutation relations for harmonic oscillators. Even the 
U(l) charge current algebra, Eq. (|2l]), is non-trivial because the right hand side, the "anomaly," has indeterminate 



sign as q • ns can be either positive or negative. Therefore, even to diagonalize Eq. (25) requires a generalized 
Bogoliubov-unitary transformation which appears not to be reducible into a product of separate Bogoliubov and 
unitary transformations. This difficulty is in contrast to that found in one spatial dimension where a simple 2 x 2 
Bogoliubov transformation is sufficient to decouple the currents associated with the left and right Fermi points! Of 
course the equation of motion for the spin currents, Eq. (|27]), is even more difficult to solve as it is non-linear. 

In the next section we will diagonalize the charge current equation of motion by resumming the perturbative 
expansion for the propagator. First we find the semiclassical solution by taking the expectation value of both sides 
of these equations and identifying 

(J(S;q)) = u(S;q) (28) 

and 

(J Q (S;q)) = 5 a (S;q) (29) 

as the amplitudes for charge and spin collective modes. Note that the collective mode amplitudes are real-valued 
in x-space because J^(S;q) = J(S; — q). The replacement of the current operators by the c-numbers u and S may 
be accomplished formally by introducing a coherent state basis that spans the space of volume-preserving geometric 
distortions of the Fermi surfaced For example, the charge collective mode coherent states arc generated by exponentials 
of the charge current operator: 

|*[u]> =exp|^ E ^' n s) ^(S;~q)}|0) (30) 

S q 

where |0) represents the quiescent Fermi liquid. A simple computation then shows that 

<*[u]|J(T;p)|*[u]> 



(*NI*N> 



u(T;p) (31) 



consistent with our definition Eq. fl28| ) . In the spin equation we also must decouple the expectation value of the product 
of two spin current operators [the third term on the right hand side of Eq. (|27j)] into the product of expectation values 
S(T;k) x S(T;q — k). This decoupling is exact in the semiclassical limit of a macroscopically occupied zero sound 
spin mode. 
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To shed light on the relationship between the semiclassical limit and the quantum regime, we present an alternative 
derivation of the collective mode equation based on the usual identification of the pole in the two-point Green's 
function, but now in the presence of a background collective mode field u(S; q). For simplicity we focus on the charge 
sector. First note that H c is diagonal in q space: the Hilbert space breaks up into a direct product of subspaces with 
different q and — q. (States with — q on the hemisphere of the Fermi surface with q • ris > 0, which we may call the 
"left" hemisphere, are coupled to states of +q on the opposite "right" hemisphere due to the indeterminate sign for 
the quantum anomaly.) Thus we may treat each (q, — q) sector separately. Now we wish to compute the retarded 
Green's function: 

r ru,i. a - (*Mk(S;g,t) J(S;-q,0)|^M) 

GWM;S;q,t) = (^frpMj m (32) 

where 

J(S; q, t) = exp[i H c t] J(S; q) exp[-i H c t] (33) 

is the current operator in the Heisenberg picture. We may choose hs ■ q > so that J(S; — q, 0) creates a particle-hole 
pair at time t = while J(S;q, t) destroys a pair at a later time t > 0. The crucial step in our calculation is to 
ignore operator ordering within the time evolution operators exp[±iH c t\. This approximation is exact so long as the 
zero mode has macroscopic occupation since in this case the errors introduced by ignoring operator ordering are small 
compared to the total energy. In other words, we should think of the operator J(S; — q) as removing just one quantum 
out of the large number of quanta that make up the macroscopic zero mode. Macroscopic occupation corresponds 
to |w(S;q)| 2 >> £1 |ng • q| which in physical terms means that there are a large number of quanta at each point in 
momentum space on the Fermi surface. (Macroscopic occupancy is possible only in the w-limit of A >> |q|; in the 
opposite q-limit the Pauli exclusion principle keeps the occupancy small.) Assuming macroscopic occupancy we have: 

exp{i H c t} J(S;q) exp{-i H c t} = J(S;p) exp{i J5[«(S;q)] t} (34) 

where E[u] is the c-number energy given by 

(*[u]|[J(S;q,t), H c ] J(S; -q, 0)|*[u]) 



£[u(S;q)] 



(*M|J(S;q,t) J(S;-q,0)|*M> 



q n s 



(35) 



Thus the Green's function is: 



G Ret ([u]; S; q, t) = |«(S; q)| 2 exp{i E[u] t} 6(t) ; (36) 



its poles in frequency space at u> = E[u] clearly correspond to solutions of the collective mode equation Eq. (|2q ) in 
the semiclassical limit. 



IV. QUANTIZED BOSONS 

In this section we calculate the Green's function in the quiescent state in the absence of macroscopic excitations. In 
this case the semiclassical approximation is inapplicable and the problem must be treated quantum mechanically. To 
simplify the calculation we restrict our attention to the case of spherical (circular in two dimensions) Fermi surfaces 
and only a single Fermi liquid parameter, the constant term Fq. Furthermore we consider only spinless fermions and 
again set the Fermi velocity equal to one. None of these simplifications is essential. 

First we write the currents in terms of boson operators that satisfy canonical commutation relations. The choice: 



J(S; q) = y/Sl |n s • q| [o(S; q) 6(h s ■ q) + a t(S; -q) 6(-n s ■ q)] (37) 

with 

[a(S;q), at(T;p)]=^ 1 C' ( 38 ) 

and 9{x) = 1 if x > and is zero otherwise, satisfies the U(l) Kac-Moody commutation relations Eq. (J2l| ) up to a 
factor of 2 which does not appear here since the fermions are spinless. The Hamiltonian Eqs. ( |l9| ) and ( |20| ) can now 
be written as H c = H + H lnt where: 
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and 



h = J2 { E e ^ s • <j) ( fis ■ q ) flt ( s ; q M s ; <*) + E e (- fis ■ <*) (- ns ■ <o flt ( s ; °( s ; -<*)} 

q ^ S S J 

fltat = E \BR E ( fis ' <0 ( fiT ' <0 V(fis-q)(fiT-q) a f (S; q) a(T; q) 

q ^ S,T 

+ 9L E (~ ns ■ q ) (~ nT ■ q ) V / (-ns-q)(-n T -q) a f (S; -q) a(T; -q) 

S,T 

+ .g ^ 0(n s • q) 9(-h T ■ q) V(n s ■ q)(-n T ■ q) a(S; q) a(T; -q) 
S,T 

+ 5 J! g (- ft s ■ q) #(n T ■ q) V(-fis-q)(fiT-q) a f (S; -q) a f (T; q) | 



(39) 



(40) 



with couplings g R = gt = g = g = /o ■ ^ will be convenient to denote a(S;q) and a(S; — q) by a#(S;q) and 

Oi(S;q), respectively the right and left moving fields. 

The generating functional for the zero temperature correlation functions is given by an integral over the coherent 
state eigenvalues Of(S; q, t) and a*(S; q, t): 



Z = J Da* T>a exp ii J dt [* a *^ ai - i? c (a*,a)] 



(41) 



where there is an implicit sum over i = L, R and the patch index S which has been suppressed. The momentum- 
frequency space propagator 

iGi(S;q,w) = <a 4 (S;q,^) oJ(S;q,£j)) (42) 
is related to the propagator of the <j) fields by: 

(&(S;q,w) &(S;-q,-w)) = - " (a<(S;q,a;) oJ(S;q,w)) . (43) 

47rns • q 

We now calculate the propagator perturbatively with the use of the bare right and left propagators: 



iG° R (S;q,u) = (oR(S;q,cj) a^(S;q,w)) 



i 



to — n s ■ q + iff sgn(w) 

and 

*G°(S;q,^) = (a L (S;q,co) a[(S;q,c)) 



(44) 



uj + ns • q + if] sgn(w) 



(45) 



The bare propagators are depicted in Fig. [y(i). 

At first order in iJ; n t there is only one connected contribution to the right two-point function and it is given by: 

iG« = (-») (»G° ) g R (n s ■ q) (<G° ) (46) 

where we have suppressed the patch, momentum, and frequency labels of the Green's function. Amputating the 
external legs, we find that the first order contribution to the self energy is just: (S; q, lu) = g R (ns • q). At higher- 
orders it is easy to see that the anomalous couplings g and ~g occur in pairs. In particular at second order there are 
two contributions to the one-particle self energy, and these are shown in Fig. [j](ii). As we build up the complete set of 
contributions to the right moving propagator, we split each scattering process, for example those shown in Fig. |^ into 
a forward scattered contribution (which involves an intermediate state in same patch S) and a remainder in which 
the boson has been scattered into a virtual state in a different patch. We can then construct the Dyson equation, 
depicted schematically in Fig. 0, where the irreducible self energy S 7 (S; q, w) comprises all amputated diagrams that 
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cannot be split into two by cutting a single bare right moving propagator. At the second order the contribution to 
the irreducible self energy is therefore: 



£( 2 >(S;q, W )=. 9 2 (n s ■ q) { £ 0( fi T ■ q) (n T • q) G° R (T; q, w) + £ 0(-n T • q) ( n_T " q) G°(T;q,c)] 



(47) 



Now we specialize to the case of two spatial dimensions. The sums over patches can be converted to integrals in the 
N — > oo limit where A — > as: 

Aj2 = k F d<f>= 2ttA(0) / dcj) (48) 
s Jo Jo 

where N(0) = lp- for spinless fermions in units where vp = 1. The second order contribution to the self energy can 
now be written more concisely as: 

S p) (S; q, w) = -/ • q) X°(q, w) (49) 

where 

x°(*) = *(o) ^ cosW ^ 

j 27r cos(</>) — x — irj sgn(cj) 
= A(0) 2 +^' 



\J x 1 — 1 yl a; 2 

= A(0) n Q (ac) (50) 

is the two-dimensional Lindhard function and x = t^t. The exact solution of the Dyson equation Fig. [3] is then given 
by: 

£ 7 (S; q, w) = J^A ( fis . q) [1 - / x (q, w )] (51) 

where 

= TT7^ T " (52) 
i + /o r'(q, w ) 

Here we see the equilibrium Fermi liquid stability criterion Fq — /gAT(O) > — 1 is necessary to keep the self energy 
non-singular in the q-limit of \x\ << 1. A little algebra then shows that the exact right moving boson propagator can 
be written in the compact form: 

i 

iG R (S ; q, w) = y- (53) 

c-n s -q|l+(|^ [l-/oX(q,w)]J 

Quasiparticle damping occurs in the q-limit when the Lindhard function has an imaginary part. In this regime we 
may write: 



Im fox(x) = Im 



= Im 



F n (x) 
i + F n ( x ) 



i - A [i-n (x)} 

A \x\ (54) 



where Aq = and the boson Green's function then reads: 



iG R (S;q,uj) = ijw - v' F h s ■ q + m s q | 



(55) 
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where the velocity is slightly renormalized from its bare value of unity: v' F — 1 + Fq{\ — Fo) . The boson lifetime 
is now finite because of scattering into different patches. Note, however, that as the self energy Eq. (|5l|) scales to zero 
as A — > it represents an irrelevant correction. In particular, the pole in the boson propagator remains unchanged 
as N — > oo. We expect this to be true generally, regardless of the shape of the Fermi surface, the details of the Fermi 
liquid parameters, or whether the fermions have spin or not. The renormalization group calculation of the second 
section therefore holds, without alteration, when the Fermi liquid interactions ZS and ZS' are turned on. 



V. FERMION QUASIPARTICLE PROPERTIES 

In the previous section we saw that Fermi liquid interactions modify the boson propagator. Though large-angle 
scattering processes were ignored, small-angle scattering processes made the boson lifetime finite. With these results 
we can use the bosonization formula Eq. (tth to infer the fermion quasi-particle lifetime. Since bosonization is carried 
out in (x, t) space we must carry out three operations. First we Fourier transform the boson propagator into real 
space. Next the exponential of the resulting expression yields the fermion propagator in real space. Finally an inverse 
Fourier transform of the fermion propagator back into momentum space allows us to extract the self energy. 

It is difficult technically to perform these steps in all generality. It will be sufficient for our purposes to first expand 
the boson propagator in powers of /o, perform the three operations on each term, and then reassemble the pieces 
to find the fermion self energy. Further, as we are interested only in the leading (second order) contribution to the 
imaginary part of the self energy, we can avoid the first of the two Fourier transforms. The real-space and time boson 
Green's function may be written as: 

iGt(S;x,t) = fa(S;x,t) ^(S;0,0) - 4> 2 {S] 0, 0)) = FT {iG^} 

= FT{ l Gf+iG^ ) +iGf + ...} (56) 

where FT represents the Fourier transform operation that converts the variables (q, uj) and to (x, t) . In the second 
line, iG<[, which is given by Eq. (fl3|), has been expanded in powers of fo. The Fourier transform of the leading 

term, FT{iG^}, is given by Eq. (|l^). Rather than Fourier transforming the first and second order corrections, we 
exponentiate this expression to obtain the fermion propagator: 

iG^S-^t) = (^(Sjx,*) ^(5; 0,0)) 



VI ^ exp{gzG,(S;x,t)} 



exp{g^ [<)+<+...]} 



(27r) 2 x ■ n s - t + iS sgn(t) 

' 1 + iFT [G$>] - i (FT [G^]) 2 + iFT [Gf] + O(/ 3 ) \ (57) 



(27r) 2 x • ns - 1 + i6 sgn(t) 

where in the last two lines we have assumed |ie_lA| << 1 and in the last line we have absorbed the factor of Air/il 2 
into Gf = {ATr/tt 2 )Gf. We are interested primarily in the imaginary part of the fermion self energy. The first order 
contribution to the boson self energy contained in iG^P is purely real and therefore does not contribute. The leading 

(2) 

contribution to the imaginary part of the fermion self energy comes from iG^' which is given by: 

iGf (S;q,u;)=-^ |lA x °( q , w ) [ w - 9|| + i v sgn(u;)}- 2 . (58) 

Since |xxA| << 1 this contribution to the boson propagator must be integrated over qj_, using Eq. (J50j) for x° an d 
we obtain: 



Urn dqxX°(S;qx,q\\,uj) = -iN(0) \u\ In j — - . (59) 

-A/2 A 



The appearance of the logarithm in this equation is peculiar to two spatial dimensions. In three dimensions the integral 
is a two-dimensional one over the two coordinates perpendicular to the Fermi surface normal and the imaginary part 
of the Fermion self energy is proportional simply to oj 2 . 
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We now take the inverse Fourier transform FT 1 . Expanding the fermion propagator as = G% + 5G^ + the 
leading imaginary contribution to the fermion propagator in (ku , uj) space is given by: 



./ 2 Aiv(o) r j , r , ,../> , y 2 - q l 



i8G^{k\\,uj)=i fC) ^ 3 I dJ \ dq\\ \uj'\ ln|- 

(60) 



(2tt)3 J_ x J_ x * A 2 

1 



(uj' — qii + irj sgn(cj')) 2 [(uj' — uj) — (q\\ — k») + irj sgn(a/ — uj)] 

The integral may be performed by complex integration; no divergences occur since all the poles in the complex qu 
plane lie either on one side of the real axis or the other unless uj' lies between and uj. Thus, except for this limited 
range of frequencies, the contour in the g|| plane may be closed without enclosing any poles. The result is: 

iSG^k^) = -I T , I 1 —7 „ a sgnH ([c 2 + (--fc||) 2 /4] In ~ *" 1 



2 (2tt) 2 [uj - k\\ +ir) sgn(w)] 2 & v ' [ L v lu ' 1 A 
+ [c 2 -(c-fc || ) 2 /4] ln^±M-|(2 W -fcu)} (61) 

and therefore the fermion self energy at this order is given by: 



lmZy>(k h oj) = - i^-^sgnH |[^ + (c-fc||) 2 /4] In i— ^ + [c 2 - ( w - k^/A] In . 



a; 



2 (2uj | . (62) 

The imaginary part of the self energy at the quasiparticle pole is the inverse of the quasiparticle lifetime. The 
location of the pole has been shifted from its bare value to uj — v' F k\\ due to renormalization of the Fermi velocity, 



l'r. 



' F = I + Fq(1 — Fq) 2^7- As a result the imaginary part of the self energy at the pole is given by: 

t v (2), 1 /o 2 ^V(0) , s / 2 , F uj2 uj 2 \ 

Im E/ ( W )U = - s g n(uj) | w In — - y j (63) 

the form of which we immediately recognize from previous work on two-dimensional Fermi liquidsEHi. This quantity 
is always negative since uj 2 « Akp. Despite the appearance of the logarithm, the weight of the quasiparticle pole, Z^, 
remains non-zero at the Fermi surface; the regular Fermi liquid interaction Fq does not destroy the Fermi liquid fixed 
point. We expect that more general regular Fermi liquid interactions, the inclusion of the spin sector, or extensions 
to non-spherical Fermi surfaces will not change this result qualitatively. 



VI. SINGULAR INTERACTIONS 

Singular Fermi liquid interactions in 2D were proposed by Anderson0 and studied perturbatively by StampBlll. 
The interaction studied couples opposite spins and diverges as k — * k': 



(k + k') ■ (k - k') 

N{oy a '- a ' |k-k'i 2 



r (kk') = — '- 0(|k| - k F ) 9(k F - |k'|) . (64) 



Note that both k and k' lie off the Fermi surface (respectively above and below it). Thus this interaction is of a more 
general sort than the type Landau originally envisaged in the phenomenological theory. The interaction diverges like 
| k jpi as k — ► k'. It is related to the regular interaction that arises at second order in a Taylor-series expansion of 

the Landau function. In three dimensions this regular interaction gives rise to a T 3 In T contribution to the specific 
heatliatt 



/(k,k') = a + fe( P -q) 2 
k- 

Ik + k'llk-k 



. (k + k') • (k - k') , 

° + ,,i +- (65) 
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for q = k — k' and p = k + k'. This interaction vanishes if k and k' lie on the Fermi surface and approach the same 
point. If, on the other hand, both momenta lie away from the surface then the interaction approaches a non-zero, but 
finite, limiting value as the two momenta converge. 

The bosonized Hamiltonians Eq. ( |l9|) and Eq. ( p2~|) generalize Fermi liquid theory in a different way: S and T lie on 
the Fermi surface but q need not be zero. Nevertheless, the interactions mentioned above havejaatural counterparts in 
the bosonized theory. The T 3 In T contribution to the specific heat is recovered in this pictureH by setting k = ks + q 
and k' = kT — q in Eq. ( |35| ) to obtain: 

V C (S, T; q) = \ IT* VF + { a + 4& + . (66) 



Note that this contribution to the specific heat is a sub-leading correction which just reflects the fact noted above that 
any q-dependence of regular Landau parameters is irrelevant to the leading-order behavior. That both generaliza- 
tions yield the same non-analytic thermodynamic behavior suggests that they are equivalent up to irrelevant terms. 
Therefore we proceed to make the same substitution in the singular interaction Eq. (|64|) of Anderson and Stamp to 
find in the spin sector: 

V S (S, T, q) - - v F S StT + — — ^ _____ (67) 

where (3 interpolates between the Anderson-Stamp interaction (/3 = 1) and the regular interaction (/3 = 2) (now for 
the spin sector). We will assume that S ^ T in the second term; otherwise the linear dispersion relation is destroyed 
at the outset and consider the scaling of the largest part of the interaction as we take N — ► oo. Observing that the 
largest contribution comes from nearest-neighbor patches S and T, and |q| < A << A, we see that the interaction 
scales as 

N 2a -P , (68) 



V N(0) 

where the exponent a was defined in section II by the equation A = kp/N a , and therefore the interaction is singular 
only if /3 < 2a < 2. 

The Anderson-Stamp interaction is singular only as S — > T. An example of an interaction which is singular for all S 
and T and which has more obvious physical significance is afforded by the Coulomb interaction. The bare interaction 
may be factorized into contributions to the three channels (ZS, ZS', and BCS). We assume that the BCS channel 
renormalizes to zero since it is repulsive. Furthermore, for small |q| the ZS' exchange channel is much smaller than 
the ZS direct channel. In this limit we find: 

7 c (S,T;q) = in- 1 « F 4 T + I 2 -^p (69) 

in three spatial dimensions with V s containing only regular interactions. Of course the identification of the bare 
Coulomb interaction with the small-angle scattering amplitude neglects the physics of screening. Nevertheless it 
would be interesting to determine the effect of the bare interaction Eq. ( |69| ) on the single quasiparticle lifetime. If 
the technical problem of performing the Fourier transforms mentioned in the previous section can be overcome, the 
bosonization method would yield non-perturbative insight into the effect of singular interactions on the Fermi liquid. 



VII. DISCUSSION 

The Coulomb interaction Eq. ( |69| ) mentioned in the previous section illustrates the difference between collective 
modes and single-particle excitations. If we neglect screening, naively substitute the Coulomb interaction Eq. (|69| ) 
into the charge collective mode Eq. (p5|), r arid compute the spectrum, we find in three dimensions a gap comparable 
to the plasma frequency; see, for exampleEj. Thus charged Fermi liquids do not support low energy collective modes 
in the charge sector. Nevertheless, we know that the single-particle spectrum remains gapless. Thus, the spectrum of 
single-particle bosonic excitations about the quiescent state, unlike the collective modes, must also remain gapless. 

The fact that the bosonized Hamiltonian separates into a sum of charge and spin parts, H = H c + H s , raises the 
specter that, as in one dimensional, the quasiparticle propagator might also exhibit spin-charge separation, even in 
the case of regular Fermi liquid interactions. Spin-charge separation in dimensions larger than one would, however, 
destroy the Green's function approach to Fermi liquid theory as the key element in that approach, the existence of a 
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pole in the single-particle Green's function with spectral weight < Z < 1, would be replaced by a branch cut and Z 
would equal zero. Fortunately this does not happen because, as we saw at the end of section IV, the location of the 
pole of the boson propagator is unchanged from its free value in the A -> limit. Consequently the spin and charge 
velocities are equal and spin charge separation does not occur. 

Finally we note that our bosonic analysis of the renormalization group flows near the Fermi liquid fixed point 
does not rely on a particular form for the fermion propagator. This is in contrast with Shankar's approach which 
assumes that the one-particle propagator always retains a Fermi liquid form; consequently non-Fermi liquid fixed 
points are ruled out from the outset of the calculation. The bosonized theory contains non-perturbative information; 
only technical difficulties prevent us from evaluating directly the non-perturbative fermion propagator. It should be 
possible to surmount these difficulties. 
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FIG. 1. Boson Green's functions, (i) Right and left moving bare boson propagators G^(S;q,a;) and G°(S; q,o>). (ii) The 
two second order contributions to the self energy which involve virtual states on the the right and left sides of the Fermi surface. 

FIG. 2. Self energy at second and third order for the boson propagator, (i) The second order contribution to the self energy 
which involves virtual states on the right side of the Fermi surface. The first diagram on the right hand side of the equation 
(with two crosses) represents scattering into and out of the same patch S. The second diagram represents scattering into and 
out of a different patch T ^ S (denoted by a dashed line with a slash), (ii) Some of the third order contributions to the self 
energy. Not shown are contributions which involve virtual states on the left side of the Fermi surface. Of the diagrams shown, 
only the first (with two dotted lines) contributes to the irreducible self energy E 7 ; the remaining three diagrams break into two 
pieces when one of the bare propagators is cut. 

FIG. 3. The Dyson equation for the self energy. The double line represents the exact one-particle boson propagator. 
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